Diffractive Excitation of Heavy Flavors: Leading Twist Mechanisms 
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Diffractive production of heavy flavors is calculated within the light-cone dipole approach. Novel 
leading twist mechanisms are proposed, which involve both short and long transverse distances 
inside the incoming hadron. Nevertheless, the diffractive cross section turns out to be sensitive 
to the primordial transverse momenta of projectile gluons, rather than to the hadronic size. Our 
calculations agree with the available data for diffractive production of charm and beauty, and with 
the observed weak variation of the diffraction-to-inclusive cross section ratios as function of the hard 
scale. 
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I. INTRODUCTION 

Diffraction is usually viewed as shadow of inelastic pro- 
cesses. This idea originated from optics helped consider- 
ably in the interpretation of data on elastic hadronic scat- 
tering. The understanding of the mechanisms of inelastic 
diffraction came with the pioneering works of Glauber [l[ , 
Feinberg and Pomeranchuk Q, Good and Walker [3[. If 
the incoming plane wave contains components interacting 
differently with the target, the outgoing wave will have 
a different composition, i.e. besides elastic scattering a 
new state will be created (e.g. see in Q). An optical 
analogy for inelastic diffraction would be a change of the 
light polarization after passing through a polarizer which 
absorbs differently two linear polarizations. 

Thus, inelastic diffraction of hadrons should also be 
treated as a shadow, it emerges due to difference of shad- 
ows produced by different inelastic processes (in this re- 
spect a study of the structure function Q of a shadow 
might look questionable). Some of the shadows from dif- 
ferent hadronic components are large, and this gives rise 
to the main bulk of soft diffractive events. However, the 
difference between shadows of two hadronic components 
which differ from each other only by either the presence 
or absence of a hard fluctuation (heavy flavors, high pt 
partons, heavy dileptons, etc.) should be vanishingly 
small. This small difference gives rise to hard diffrac- 
tion. It can be as small as 1/Q 2 , where Q is the hard 
scale, and then the amplitude squared leads to higher 
twist terms in the diffractive cross section. However, in 
a non-abelian theory like QCD the difference between 
shadow amplitudes may be larger, ~ resulting in a 
leading twist behavior. 

A proper example is the leading twist diffractive Drell- 
Yan reaction where simultaneously large and small 
size projectile fluctuations are at work. Due to general 
properties of diffraction a single quark cannot radiate 
diffractively a photon or a dilepton (or any colorless and 
point-like particle) in forward direction (when the recoil 
target proton has transverse momentum pt = 0) Q • Let 
us consider two Fock components of a physical quark, a 



bare quark \q), and a quark accompanied by a photon, 
|<77*). In both Fock states only the quark is able to inter- 
act. Therefore, the forward diffractive amplitude given 
by the difference of these two elastic amplitude integrated 
over impact parameter, vanishes. 

Nevertheless, a qq dipole (or any hadron) can radiate 
diffractively in forward direction. Indeed, in this case 
the two Fock states \qq) and |<7<?7*) have different cross 
sections and the difference is, 

A dlff (qq +p -> qq~f* + p)\ PT=a cx <Jq q (R) - <J qq (R + r) 

cx f-R + 0{r 2 ). (1) 

Here R is the qq transverse separation; r ~ 1/Q is the 
small shift of impact parameter of the quark radiating the 
heavy photon; a qq {R) is the dipole-proton cross section 
which we assume for simplicity to be quadratic in R. It 
is demonstrated in [|| that integration over azimuthal 
angle does not terminate the leading term (r • R) due to 
the convolution with the final state wave function. Thus, 
the amplitude has a leading twist scale dependence, 

The origin of leading twist behavior of heavy flavor 
diffractive production has some similarities, but also dif- 
ferences. In this paper we classify different mechanisms 
of diffractive excitation of heavy flavors, with the pur- 
pose of identifying leading twist terms. It is instructive 
to compare them with the well studied case of diffractive 
deep- inelastic scattering (DIS). 



A. Diffractive production of heavy quarks in DIS 

The fraction of DIS events with diffractive excitation 
of the virtual photon has been found in experiments at 
HERA to be nearly scale independent This might 
have been a surprise, since diffraction, which is a large ra- 
pidity gap (LRG) process, is associated with a Pomeron, 
i.e. at least two gluon exchange. If each of the gluons 
has to resolve a hard photon fluctuation, 7* — > qq, of a 
size 1/Q 2 (Q is the photon virtuality), then to do it twice 
costs more and should lead to a cross section as small as 
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l/Q 4 . However, such an expectation contradicts data. 
It turns out that in this case the main contribution for 
transversely polarized photons comes from rare fluctu- 
ations of the photon. These fluctuations correspond to 
aligned jet configurations @ which happen rarely, but 
are soft and interact strongly [Ic| |. 

Although the cross section of diffractive dissociation 
7* — > qq seems to behave as leading twist, a s d("f*P — > 
qqp ~ l/Q 2 ), in fact this is a higher twist process. In- 
deed, if we impose the hard scale to be the heavy quark 
mass rug ^> Q 2 , the same process behaves as one would 
expect for a higher twist, 

<r sd (i* QQ + p) oc l/m%. (2) 

The real leading twist behavior emerges from more 
complicated photon fluctuations which contain at least 
one gluon besides the QQ, 7* — > QQG (e.g. see in 
@,E|)- Such fluctuations are characterized by two sizes, 
one controlled by the hard scale, which is small, either 
r ~ l/Q 2 or r ~ 1/toq. Another size, the mean quark- 
gluon separation, is large and depends only logarithmi- 
cally on the scale. This fact gives rise to the leading twist 
behavior of diffraction: one of the ^-channel gluons has 
to resolve the small (l/Q 2 or l/m|>) size, while another 
gluon may interact with the large QQ — G dipole. This 
corresponds to diffractive excitation, which is different 
from ©, 

<Jsd(l* +p —> QQg +P) oc l/mg. (3) 

Gluon radiation in the final state is here essential for 
the leading twist behavior. Indeed, although higher Fock 
components, like QQg, also contribute to ((2|), but this 
does not not change its higher twist scale dependence. 

Notice that diffraction is closely related to nuclear 
shadowing, since both emerge from the unitarity rela- 
tion as a shadow of inelastic processes. A direct relation 
between diffractive excitations of the beam and nuclear 
shadowing was first found in [12] , and is known as Gribov 
inelastic shadowing. In our case, the leading and higher 
twist contributions to shadowing are related to the same 
types of fluctuations of the photon, QQG... and QQ re- 
spectively (HQ. 



B. Diffractive hadroproduction of heavy quarks 

One may treat diffraction in DIS as a way to measure 
the partonic structure of the Pomeron Q. Having this 
kind of information one may try to predict other hard 
diffractive processes assuming factorization. However, 
attempts to apply QCD factorization to hard diffrac- 
tion failed by an order of magnitude when one compares 
diffraction in DIS and in hadronic collisions [151 ]. 

There are many reasons for this breakdown of factor- 
ization. The first one is pretty obvious, and comes from 
the absorptive or unitarity corrections, which has been 



known since the era of Regge phenomenology. These ef- 
fects cause the suppression of any LRG process, except 
elastic scattering. In the limit of unitarity saturation 
(black disk) the absorptive corrections may completely 
terminate the LRG process. Actually, this almost hap- 
pens in (anti)proton-proton collisions, where unitarity is 
nearly saturated at small impact parameters (l6| . The 
suppression factor, which is also called survival proba- 
bility, changes the diffractive cross section by an order 
of magnitude. Although hard reactions hardly make 
any shadow, the strength of absorptive corrections in 
hadronic collisions is controlled by the soft spectator par- 
tons (see Sect. I VII which are absent in the case of diffrac- 
tion in DIS. This is why factorization is severely broken. 

Another source of factorization breaking is the differ- 
ence between the mechanisms of diffractive QQ pair pro- 
duction in DIS and in hadronic collisions. In both cases 
the Pomeron (i.e. two or more gluons) can be attached 
directly to the produced heavy quarks, and this part of 
the interaction is subject to factorization. In a hadronic 
collision, however, the Pomeron can be attached simulta- 
neously to the projectile gluon and to the heavy quarks. 
In other words, the heavy pair, which has a lifetime sub- 
stantially shorter than the projectile gluon in the incom- 
ing hadron, may be produced during the interaction pjj • 
This part, called coherent diffraction [l8j], causes another 
deviation from factorization. 

This part of diffraction was modeled in 17] by a col- 
orless projectile gluon diffractively dissociating into QQ. 
The cross section was found to be leading twist, oc l/mg. 
The confusion caused by the graphic presentation (Figs. 3 
and 5 of [17| ) motivated a more complete calculation in 
fl9| of the full set of Feynman graphs corresponding to 
diffractive dissociation of a colored gluon, g+p — > QQ+p. 
These authors found this process to be a higher twist, like 
in photoproduction Eq. Indeed, since in the dissoci- 
ation g — > QQ all transverse distances between the qluon 
and quarks are of the order of I/vtiq, the cross section 
of this process mediated by Pomeron exchange must be 
oc l/nig. Then this part of diffraction violating factor- 
ization has the same toq dependence as the factorized 
one, i.e. both are higher twists. 

A new source for breakdown of factorization was found 
in [6( for diffractive Drell-Yan processes. It turns out 
that in this case the participation of soft spectator par- 
tons in the interaction with the Pomeron is crucial and 
results in a leading twist effect. A similar mechanism for 
hadroproduction of heavy quarks is under consideration 
in the present paper. It is related to the processes, 

c s d(g+p— > qQQ + p) oc l/m 2 Q . (4) 
°sd{g + p-> gQQ + p) oc l/m 2 Q . (5) 

Just as in leading twist diffraction in DIS, Eq. ([3]), these 
processes are associated with two characteristic trans- 
verse separations, a small one, ~ I/tuq, between the Q 
and Q, and a large one, either ~ between q and 

QQ in (|1]), or ~ l/«i g between g and QQ in ([5]). Here 
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m q ~ Aqc d and the several times larger m g (see [lQ, [20| 
are the effective cut-off parameters which take care of 
the nonperturbative interactions of quarks and gluons 
respectively. 

Somewhat similar, but nevertheless different observa- 
tions were made in Ref. [Hj]. Namely, similar to Drell- 
Yan diffraction Q the large hadronic size enters the 
leading twist diffractive amplitude due to the interac- 
tion of the projectile remnants with the target. There 
it was concluded that theoretical predictions cannot be 
certain since we are lacking reliable information about 
the hadronic wave function. However, the main leading 
twist contributions Eqs. Q-© under consideration in 
the present paper are independent of the structure of the 
incoming hadron. They correspond to diffractive excita- 
tion of an individual parton via the so called production 
mechanism (see Sect. ITT]) . 

Interaction with spectators is also known [12] to cause 
considerable effects in the azimuthal single-spin asymme- 
try in semi-inclusive pion leptoproduction. In this case 
the outgoing quark experiences final state interactions 
with remnants of the proton. 



C. Intrinsic heavy flavors in the proton 

Production of heavy flavors at large Feynman x f has 
been always a controversial issue, even in the simple case 
of inclusive processes. In the perturbative QCD approach 
based on QCD factorization, inclusive heavy quark pro- 
duction is described as glue-glue fusion, and the rapidity 
distribution of produced heavy flavored hadrons is con- 
trolled by the gluon distribution in the colliding hadrons. 
The gluon density steeply vanishes towards Xf = 1, ap- 
proximately as (1 — xp) 5 '■ Convolution with the frag- 
mentation function (assuming factorization) makes this 
behavior at xf — > 1 even steeper. On the other hand, 
the end-point behavior is dictated by the general result 
of Regge asymptotics, 



da 



dxF dt 



OC (1 - SBf) 1-20 ***) 



(6) 



xp — >1 



where an(t) is the Regge trajectory corresponding to the 
i-channel exchange of a heavy flavored meson or baryon, 
which depends on the quantum numbers of the projec- 
tile and the produced heavy flavored hadron. Apparently 
this has little to do with the gluon distribution function. 
The same problem appears in the Drell-Yan reaction at 
xf —* 1) as is seen in data (23[. Therefore, one should not 
rely on QCD factorization at large xp — ► 1. In fact, in 
this kinematic region several mechanisms breaking fac- 
torization are known [24j |. 

An excess of heavy flavored hadrons at large xp is 
expected to be an evidence for the presence of intrin- 
sic heavy flavors in the projectile hadron [2^, [2(|. No- 
tice, however, that any mechanism must comply with the 
Regge behavior of Eq. © independently of the details of 



hadronic structure, i.e. must be the same with or with- 
out the presence of intrinsic heavy flavors. Besides, to be 
confident that an excess of heavy flavor is observed, one 
must be able to provide a reliable theoretical prediction 
for the conventional mechanisms at large Xf- This is a 
difficult task in the situation when QCD factorization is 
broken. 

The observation of diffractive production of heavy 
quarks may provide a better evidence for intrinsic heavy 
flavors. This is expected to be seen as an excess of 
diffractive production compared to the conventional ex- 
pectation. The latter, therefore, must be reliably known. 
However, an observed signal might be misinterpreted if 
it is compared with calculations assuming that factoriza- 
tion holds for hard diffraction 0, [27J ■ This is not cor- 
rect as was explained above and is confirmed by direct 
calculations shown below. A good signature for the con- 
tribution of intrinsic heavy flavors would be the sharing 
of longitudinal momentum in the diffractive excitation. 
The intrinsic heavy quarks should carry the main part of 
the momentum, which would be unusual for the conven- 
tional mechanisms (see Sect. lVII)) . However, this requires 
the observation of both heavy quarks, which is difficult. 

Notice that diffractive production of heavy flavors also 
creates a large background for searches of Higgs bosons, 
which can be produced at large xf from intrinsic heavy 
flavors 28]. Moreover, one also needs to know the rate 
of direct production of heavy flavors while searching for 
physics beyond the standard model. 



D. Outline of the paper 

We present below a calculation for diffractive produc- 
tion of heavy quarks in proton-proton collisions. For this 
purpose we rely on the dipole approach, which is an al- 
ternative phenomenology to the parton model formalism 
based on QCD factorization. The key ingredient is the 
universal dipole cross section introduced in [29| . which 
is the total cross section, aq q {rT,x), of interaction of a 
quark-antiquark dipole of transverse separation tt with 
a proton. The Bjorken variable x depends on rx and 
the dipole energy. Notice that this is essentially a target 
rest frame description, interpreting the beam hadron as 
a composition of different Fock states. Those light-cone 
fluctuations are assumed to be " frozen" by Lorentz time 
dilation during the interaction. Thus, the cross section 
is a sum of the dipole cross sections for different Fock 
states. At the same time, these dipole cross sections are 
averaged over the properties of the target. 

The paper is organized as follows. In Sect. [IT] we start 
with inclusive production of heavy quarks and classify 
the different mechanisms. The first one, which we call 
Bremsstrahlung, is similar (except for the couplings and 
color factors) to the Drell-Yan mechanism. It arises from 
the interaction of the radiating color charge, similar to 
bremsstrahlung in QED. The second one, called Produc- 
tion mechanism, is related to direct interaction of the 
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radiated virtual gluon or the heavy quark pair with the 
target. The important observation of this section is the 
smallness of the interference between the two mecha- 
nisms. 

In Sect. IIIII we calculate the forward amplitudes for 
diffractive production of a QQ pair in quark-proton col- 
lisions. The two mechanisms, bremsstrahlung and pro- 
duction, are found to have different scale dependence. 
The former is a higher twist effect, while the latter is 
leading twist and dominates the diffractive cross section. 
We generalize these results to diffractive gluon-proton 
collisions in Sect. (|IVp . 

The next step is the calculation of diffraction in proton- 
proton collisions. In this case the Pomeron, which is a 
multigluon exchange, can interact with active (radiat- 
ing) and spectator partons simultaneously. This possi- 
bility gives rise also to a leading twist contribution for 
the bremsstrahlung mechanism. 

In sect. [Vl]we estimate the suppressing effects of uni- 
tarity saturation, which is present in pp collisions at high 
energy. Our results are close to other estimates available 
in the literature. The suppression caused by absorption 
ranges from an order of magnitude at the Tevatron down 
to a few percent at LHC energies. 

Numerical calculations are performed in Sect. IVII1 re- 
lying on the phenomenological dipole cross section, well 
fitted to HERA data on the proton structure function. 
We present the results for energy dependent diffractive 
cross sections of different heavy flavor production, x\ 
distribution (fractional momentum carried by the heavy 
quark), and transverse momentum dependence at differ- 
ent energies. We found that the bremsstrahlung mech- 
anism, although a leading twist, gives a negligibly small 
contribution. The main part of the cross section cor- 
responds to the production mechanism in collisions of 
projectile quarks and gluons with the target proton. 

In Sect. lVIlH we review the available, rather scarce data 
on diffractive production of heavy flavored hadrons. The 
data agree quite well with our calculations. Our results 
and observations are summarized in Sect. HXl 



II. 



INCLUSIVE PRODUCTION AT FORWARD 
RAPIDITIES 



A comprehensive study of inclusive heavy flavor pro- 
duction within the light-cone dipole approach was per- 
formed in [3(3] (see also review [3l|). However, those 
calculations were focused on heavy flavors produced at 
mid rapidities, where gluon-gluon fusion is the dominant 
mechanism. Here we are interested in heavy quarks pro- 
duced in the projectile fragmentation region, and our ul- 
timate goal is diffractive production. In this case the 
interaction with valence quarks must be included, and 
this needs more elaborated calculations. 

First, we calculate the amplitude of inclusive cross sec- 
tion for the production of a heavy pair QQ in a quark- 
proton collision, q + p — ► qQQX, in one gluon approxi- 



mation, as is shown in Fig. [T] 



q(p,) 




FIG. 1: One gluon graph for heavy quark pair production by 
a quark in inelastic collision 

In what follows we use the following notation (see 
Fig. [1]): pi 2 are the 4-momenta of the projectile quark 
q and the ejected quark q' respectively; fci t 2 are the 4- 
momenta of the produced heavy quarks Q and Q re- 
spectively. For momenta combinations: k = k\ + k2] 
a = k+/p+; 13 = kf/k+- M 2 = (h + fc 2 ) 2 ; Q 2 = -{pi - 
P2) 2 > 0; 7? = ap 2 — (I — a)k is the relative transverse 
momentum between q' and (QQ); K = (1 — 0)ki — (3k 2 
is the relative transverse momentum between the heavy 
quarks. 

To make further progress we switch from transverse 
momenta to impact parameters, which require further 
definitions: r, ri, r*2 are the transverse separations 
within the Q — Q, q — Q and q — Q pairs respectively; 
p = r — 0f\ — (1 — /3)r*2 is the distance between q' and 
the center of gravity of the QQ pair; s = ri — fa is the 
QQ transverse separation. 

Since we are interested in QQ production in the projec- 
tile fragmentation region, the rapidity interval between 
the q' and QQ is assumed to be short. At the same time 
the rapidity interval between the QQ and the target is 
long (at high energies) and is filled by radiated gluons. 
The bottom blob in the graph in Fig. Q] includes all those 
gluons, while the upper part of the graph can be calcu- 
lated using the Born approximation which is rather ac- 
curate for the projectile fragmentation region. This part 
is represented by the Born graphs depicted in Fig. [51 

The sum of this five amplitudes can be spiffed into 
two classes which we assign to bremsstrahlung (Br) and 
production (Pr) mechanisms, as is described in detail in 
|Appendix A| 



i=l 



M 



T(L) 



M 



T(L) 



(7) 



where the upper indexes T and L correspond to trans- 
verse or longitudinal polarization of the virtual gluon ra- 
diated by the projectile quark. 

The transverse bremsstrahlung amplitude, Mg r , cor- 
responds to the following combination of graphs depicted 
in Fig. [2] (see I Appendix A|I), 



M T Br = Ml 



Mi 



Q 2 



M 2 + Q 



2 M i 



(8) 
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Since the effective dipole cross section is a func- 

tion of one variable, for the sake of convenience we use 
here the light-cone distribution amplitudes in mixed im- 
pact parameter - transverse momentum representation, 

$Br(fla)=5Z*l(#*2(2) + *l(d*2(2) , (13) 



where the sum is over the transverse polarizations A of 
the radiated virtual gluon; the amplitude of longitudi- 
nal gluon radiation is labeled by the index L. Further 
notation is, 



V^UTp) t 



2tt 



|z(2 - a)e x ■ V + a [a x e A ] • V 



ia 2 m q [a x e\] ■ n) Xin K q (tp) 



(14) 



FIG. 2: Feynman graphs contributing to inclusive production 
of a heavy quark pair. 



It describes bremsstrahlung of a transversely polarized 
heavy gluon which dissociates into QQ. The remain- 
ing transverse amplitudes are combined into the second 
group, the amplitude of QQ production via direct inter- 
action with the heavy quark pair, 



Mi 



M 



AP + Q 



2 M i 



Mj + Mi 



(9) 



The procedure of grouping the longitudinal amplitudes 
is more complicated and is described i nJAppendix A| 2. 
It leads also to the structure Eq. \f§, where the 
bremsstrahlung and production longitudinal amplitudes 
are defined in (1A.21|) . 

Thus, the inclusive cross section has the following 
structure, 



da(qp -> QQX) da 1 



da 1 



where 

da T ( L Hqp 



da d(3 



QQX) _ da 



da dp da dp 



(10) 



T(L) 
Br 



da 



T(L) 
Pr 



da 



T(L) 
Int 



da dp 



da dp da dp da dp 



The last term here corresponds to the interference of the 
amplitudes Mg r and Mp r . 



A. Bremsstrahlung mechanism 

The first term in the cross section Eq. pd|) can be 
calculated in the framework of the light-cone dipole ap- 
proach as follows, 



da Br (qp ^ QQX) 
da dp 



(2tt): 



■d 2 P \<5>Br{p,K)Y 



Si(p); 
(12) 



2tt 



X \ 2(1 - a)M Xin K (rp) , (15) 



which are the light-cone distribution amplitudes for the 
quark-gluon Fock state with transversely polarized gluon 
(polarization A) and longitudinal gluon respectively. The 
spinors Xin and Xf correspond to the initial and final light 
quarks. The running QCD coupling a s (l/p) is taken at 
virtuality ~ 1/p 2 - 

The distribution amplitudes for a QQ fluctuations of 
transversely and longitudinally polarized gluons respec- 
tively, read, 



*a(«) = (t>' i {m q v-ex + (l-2p)(a-n)(ex-K) 

a s (n) 



+ i [ex x n • k 



(16) 



%%{k) = 0t 2/3(1 - P)Ma -ri(j) 



a s {n) 

n Q + K 



(17) 



where 
tively: 



are the spinors of the heavy Q and Q respec- 



t 2 = (1 - a)M 2 + a 2 m 2 ; 



M 



m Q + k2 
P(l ~ P) 



m q and uiq are the light and heavy quark masses re- 
spectively; ex is the polarization vector of the transverse 
gluon; n is the unit vector aligned along the direction of 
the projectile quark; a are the Pauli matrices. 

The dipole cross section S 1 ( p) in (fl2l corresponds to 
gluon radiation by a quark [7J, |32|, |33| . It has the form 
of the cross section for a gluon-quark-antiquark color- 
less system with transverse separations pfor gluon-quark, 
(1 — a)p for gluon- antiquark, and ap for quark- antiquark, 



= | W(P) F a [(1 a)p\} - I a(ap) 



(18) 



where a(p) is the cross section of interaction of a qq dipole 
with transverse separation p on a proton. 
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B. Production mechanism 



C. Bremsstrahlung-Production interference 



The second term in (|lip is represented in a form similar 
to d, 



■d 2 s\<t> Pr {s,p 2 )\ 2 Z 2 {s), 



d(j Pr (qp — ► QQX) _ f d 2 p 2 



da df3 



. (19) 

with distribution amplitudes having structures similar to 

^Pr{s,p 2 ) = J! *3(S)**(&) + *3 (3*4 (ft) , (20) 



where 



= ^{m Q c^-i(l-2/?)(c^ri)(e s -V) 
+ [e x xn]-v| 0X O («) ; (21) 



= ^ 2/3(1 -P)Q4>K {es) 



(22) 



are the light-cone QQ distribution amplitudes in impact 
parameter representation. In the above expressions, 



e 2 = (3(1 ~ P)Q 2 + m% ; 
2 _ pi + a 2 m 2 



(23) 



Notice that differently from the bremsstrahlung ampli- 
tudes, here the gluon virtuality Q is not equal to the 
effective mass of the QQ. 

Correspondingly, the transverse and longitudinal dis- 
tribution amplitudes for quark-gluon fluctuation in mo- 
mentum representation read, 



= X/{( 2 " a ) (^A -V2) + ia[a x e x ] ■ p 2 

1 



-I- ia 2 TO 2 [<? x ej] • n} Xi 



9 9 i 2 ' 



(24) 



*i(p2) = 2(1 - a) / 2 
a z m^ + p^ 



(25) 



The dipole cross section £2(5) m ([19)) corresponds to 
gluon decay into a quark-antiquark pair [33, 132I ] . Similar 
to (fl8|) it also has the form of a cross section for a gluon- 
quark-antiquark colorless system, but with transverse 
separations s for the quark-antiquark, (3p and (1 — [3)p 
for gluon-quark and gluon-antiquark, respectively. 

Sa(S) - I + a [(1 - /?)*]} - 1 tr(S) . (26) 



The third interference term in (jlip reads, 



da d/3 



(27) 

In this case the effective dipole cross section is function 
of two coordinates, 

£i„t(p,s) = l{a(p,-p£)+a[(l-a)p,(l-p)g\ 

- a[p,(l-/3)s]}-±{a[p,(l-0)g\ 

+ a [(1 - a)p, -/3s] - a(p, -/3s)} , (28) 



where 



a(p, s) = a{p + s) - a(p) - cr(s) 



(29) 



For this reason we use in (|2T|) the distribution functions 
$ fully in impact parameter representation. They are 
related to the mixed representation functions given in 
Eqs. (T3J) and $ZU§ as, 



4> 



1 



t(l; 

Pr 



(2tt) 2 

1 

(2^ 



d 2 Ke 



d p 2 e" 



/_,, (p,«0 ; (30) 

Pr (P,P2) (31) 



This interference term, Eq. (|27p . turns out to be sup- 
pressed in the cross section compared to the first two 
terms in (jlip . Indeed, performing the integrations in 
Eqs. (Tj"2")) and Q19p . one arrives at cross sections which 
expose the leading twist scale dependence both for the 
bremsstrahlung and production mechanisms, 



d<if 



do 



Pr 



1 



da d/3 da d(3 m 2 -, 



(32) 



Indeed, the mean transverse separations squared, which 
are controlled by the distribution amplitudes Eqs. (|14p . 
(fT5| and Eqs. fl2TJ), are (p 2 ) ~ 1/e 2 ~ l/m 2 Q and 

(s 2 ) ~ 1/t 2 ~ l/m/Q, respectively. At such small sep- 
arations the cross sections Eqs. (|18l) and (|26p behave in 
accordance with color transparency [29j, oc p 2 and 

E 2 (s) oc s 2 , confirming (p{2")) . This result could be ex- 
pected, since it also follows from the QCD factorization 
scheme. 

As for the interference term, it follows from Eqs. (|T3|) 
and (|20p that the product of distribution amplitudes 

in (|27|) . ($> Pr (p, s) §Br{p, s)^, contains only coordinates 

squared, p 2 and s 2 , but not their product, p ■ s . At the 
same time, since the dipole cross section at small sepa- 
ration behaves as <r(r) = C r 2 , the combination of cross 
sections, Eq. (|29|) . reads, 



a(p, s) = 2C p ■ s . 



(33) 
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Thus, the interference term Eq. ([25]) vanishes after inte- 
gration over azimuthal angle. 

In a more realistic model with a dipole cross section 
leveling off at large separations, Eq. (f7"Tj) . the result is 
not zero, but is a higher twist effect, 



da d/3 



1 



rn 



%Rl(x) 



(34) 



where Rq(x) is defined in (|71[) . Thus, the interference 
term in the cross section of inclusive heavy flavor pro- 
duction is vanishingly small. This observation turns out 
to be valid for diffractive production as well. Therefore, 
we safely neglect the interference in what follows. 



III. DIFFRACTIVE PRODUCTION OF HEAVY 
FLAVORS IN QUARK-PROTON COLLISIONS 

Here we calculate the cross section of diffractive exci- 
tation of a projectile quark resulting in the production of 
a heavy quark pair, 

q + p — > qQQ + p , (35) 

as is shown in Fig. [3] This picture represents numerous 




FIG. 3: Diffractive production of a heavy quark pair in a 
quark-proton collision 

Feynman graphs which are equivalent to a rather sim- 
ple factorized form of the diffractive amplitude in the 
light-cone dipole representation. This formalism was de- 
veloped for diffractive gluon bremsstrahlung in 0]. 

Similar to ©-([5]), the diffractive amplitude A qp can be 
splitcd into two parts, 



A«* = A q £ r 



(36) 



where the bremsstrahlung and production amplitudes are 
defined as follows. 



A. Bremsstrahlung mechanism of diffraction 



The bremsstrahlung amplitude in (|36|) has the form, 



The distribution amplitude $g r (p, s) in coordinate space 
is given by Eq. [30|) . 



*Br(£s) = ^ L Br {p,s) + ^>l r {p,s) 



^a s (l/s)a s {\/ p) 



(2tt)2 

x [4(1 - a)^-Ki{m Q x)xfXin4>^n(j) 

+ a x(p) &a(V s ) \ K (m Q 3 
x p 



(38) 



where 



o-xip) = X/{*( 2 - a)ex ■ p + a[a x e A ] • pjxin] (39) 



bx(\7 s ) 



4>^m q a ■ e\ + [e\ x ft] ■ V s 
i(l-2/3)(o ! -n)(e A -V s )}0; 
(1 - a)p 2 




(3(1 - p) 



(40) 
(41) 



The imaginary part of the amplitude of diffractive 
production of heavy flavors is calculated employing the 
generalized optical theorem (Cutkosky cutting rules) 
|34| . It is illustrated in Fig. [4] for the example of the 
bremsstrahlung mechanism, and is used in the same 
way for other mechanisms in what follows. The ampli- 



Im 




+ 



+ 



+ 



TV 



(37) 



FIG. 4: Cutkosky rules for the imaginary part of the ampli- 
tude of diffractive QQ production in quark-proton collision. 
Only bremsstrahlung mechanism is considered for this exam- 
ple. 

tudes shown in the graphs on both sides of the unitarity 
cut (vertical dashed line) are on-mass-shell physical pro- 
cesses, and the final state excitations of the proton are 
summed up employing completeness. After integration 
over impact parameter b, this procedure leads to the fol- 
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lowing form for the effective cross section Ei in ([57)) , 



S) - Q + ^ 7*r«) {a [p - (1 - 0)5] 
- a [p + 0g\ - a [(1 - a)p - (f - 0)5] 
+ a [(1 - a)p + + \ r q a r?{a[p\ a[(l - a)p\ 

+ CT [(l-/?)3-i CT [(l- a )p- (1-/3)*] 



1 



a[p-(l-/3),s]-a 



+ -a[(l-a)p-/3sl + -a[p + /3s]} 



(42) 



Here t% — \\ q a and = ^A^, where A* and \® are the 
Gell-Mann matrices acting on the color spaces of light 
and heavy quarks respectively. 

The corresponding bremsstrahlung term in the diffrac- 
tive cross section reads, 



S 2 (P, t) = + - r«r« J {a [p + 05] 
- o- [p - (1 - /3)s\ + a [(1 - a)p + /3s\ 

- a [(1 - a)p - (1 - /3)5]} + ^ r^a[s\ + 

+ I r^{a[s\ - a[(l /3)5] - a[(3g\ + a(p) 

- a[{l-a)p\+ l -a[{l-a)p-{l-(3)s\ 
+ l a [(l-a)p + ^-^a[p- (1-/9)5] 



-a[p + /35]} 



(46) 



Eventually we arrive at the production cross section, 



da Pr (qp -> QQgp) 



(it' <ia d/3 



t'=o 



256tt 



d 2 pd 2 s 



|$ Pr (p,s)| Tr 9 Tr Q E 2 (p, s)E^(p, s) 



(47) 



da Br (qp -> QQqp) 



dt' da dp 



t'=o 



256tt 



d 2 pd 2 s 



x |i> Sr (p,s)| 2 Tr,Tr Q E^^Q?, s) , (43) 



where the traces Tr g and Ttq are taken over the Gell- 
Mann matrices corresponding to light and heavy quarks 
respectively. 



B. Production mechanism of diffraction 



The production term in (|36p reads, 



Af r = ^J d 2 pd 2 se^ + ^ s <S> Pr {p,s)t 2 {p,s) , (44) 



where 



y/a s (l/ S )a s (l/p) 
(2tt) 2 



$Pr(p, S) = ^ r (p, S) + $£ r (p, S) 

|4(1 - a) T ^-K 1 (m Q x)xfXin(t> * vncj> 

*\{p) b x {V s ) \ [K a (m Q s) - K (m Q x)} , (45) 



The functions a\(p), b\(W s ) and x arc defined in (|3"9")) - 



C. Scale dependence 

In the bremsstrahlung distribution amplitude, 
Eq. (j38j) . both mean separations are controlled by the 
hard scale, 



(P 2 



1 



(48) 



Therefore, the corresponding term Eq. (|43|) in the diffrac- 
tive cross section is as small as 1/toq, i.e. it is a higher 
twist effect. 

On the contrary, in the production mechanism only the 
QQ separation is small, (s 2 ) ~ V m Q- The mean sepa- 
ration between the light quark q' and the QQ according 
to ([24|) - ((25l) is large, (p 2 ) ~ 1/m 2 . However, the effective 
cross section Eq. (|46| cannot be large. Indeed, at small 
s — > Q it vanishes as 

E 2 (p,s)cxs-p . (49) 

This result is similar to what was found in [6] for the 
diffractive Drell-Yan reaction, which also probes simulta- 
neously large and small distances. This is very nontrivial, 
since in the case of the Drell-Yan reaction that property 
is due to the Abelian nature of the radiated particle. Now 
we have a non- Abelian radiation, but arrived at the same 
feature. 

It is interesting to notice that while the for- 
ward Abelian radiation by a quark is forbidden, the 
bremsstrahlung part of the diffractive radiation of a QQ 
pair, although is not zero, turns out to be quite sup- 
pressed. 

For further calculations we can safely neglect the 
higher twist contribution of the bremsstrahlung mecha- 
nism to the quark-proton diffractive amplitude, and keep 
only the leading twist production term, Eq. (|4"7| . 



9 



IV. DIFFRACTIVE PRODUCTION OF HEAVY 
FLAVORS BY A GLUON 



Not only quarks, but projectile gluons can also be 
diffractively excited producing a heavy QQ pair, as is 
illustrated in Fig. [5j This contribution is important if 




FIG. 5: Diffractive production of a heavy quark pair in a 
gluon-proton collision 



the QQ pair carries a small fraction of the initial beam 
momentum. 

The calculations are similar to those that we per- 
formed in Sect. IIIIBI Here we neglect the higher twist 
bremsstrahlung contribution. The leading twist ampli- 
tude of diffractive QQ pair production by a projectile 
gluon, g + p — > QQX + p, reads, 



The effective cross section in (f5"U| is given by 

S 3 (p, S) = (y-^ab + ^dabcT^j ~ Ol)p- (1 - f3)s] 

+ a[p-(l-0)s]-a[(l-a)p + /3s] 

- (r[p + Ps]} + jf abd T d {2a[(l-a)p\ 

- a[(l- a)p + [3s] - a[(l - a)p- (1 - p)3\ 

- 2a(p)+a[p + [3s]+a[p~(l-P)s] 

- 2a(f3s)-2a[(l-P)s]+-a(s)}. (53) 



The indexes a, b correspond to the initial and final state 
polarizations of the gluons, and d a b c and f a b c are the 
structure constants. 

Eventually, the differential cross section of diffractive 
production of a QQ pair in gluon-proton collision by the 
production mechanism (leading twist) turns out to be 



da Pr (Gp -> GQQp) 



dt' da dp 



256?r 



d 2 pd 2 s 



<S>f p) (p, s) Ttq ( E 3 (p, mUp, a) 



(54) 



Here the brackets (...) indicate the averaging over initial 
and summing over final color indexes of the gluons. 



Af r p) = ^J d 2 pd\se^^^\p^Mp^ , 

(50) 

where 

®p r (P>s)= ^2 \A^-Ki{m Q x)e f ■ e in 

+^ a x{p)b\(V s )^[K n (m Q s) - K (m Q x)]} 

(51) 

Here 

a x {p) = 2i[(l - a)(e f ■ e in )(ex ■ p) + a(e in ■ e\)(ef ■ p) 
- a(e f ■ e x )(ei ■ pj] , (52) 

b\(V s ) is given by Eq. (|40|) . x is given by (|4T|) . e m and e/ 
are the polarization vectors of gluons in the initial and 
final states respectively, and e\ is the transverse polar- 
ization vector of the virtual ^-channel gluon. The first 
and the second terms in the r.h.s. of Eq. (|5ip correspond 
to exchange in the t-channcl of longitudinally and trans- 
versely polarized gluons respectively. 



V. DIFFRACTIVE PROTON-PROTON 
COLLISIONS 

Now we consider a large rapidity gap single diffractive 
process, 

p + p^ QQX +p, (55) 

where one colliding proton remains intact, and the debris 
of the other proton contains a heavy QQ pair. Similar to 
the quark-proton collision, the amplitude of this reaction 
can be splited into bremsstrahlung and production parts, 

App = A™ + A p P p r , (56) 

which arc described below. 



A. Leading twist bremsstrahlung contribution 

The first bremsstrahlung term in (|56|) is represented 
by graphs depicted in Fig. [SI We neglect the higher twist 
term corresponding to attachment of both t-channel glu- 
ons to the same valence quark in the projectile proton. 
More diagrams are generated by permutation of the pro- 
jectile quarks. 
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FIG. 6: Diffractive production of a heavy quark pair in a 
proton-proton collision. More graphs emerge after quark per- 
mutation. 



The diffractive amplitude corresponding to the graphs 
in Fig. [5] has the form, 

A%. = [ d 2 pd 2 sT\d 2 r l dx l e^ + ^ 

16tt J xx 

X *J,(fi,f 2 , f 3 ,Xx,X 2 ,X 3 ) Xg^p, S, fx, f 2,1%) 

x ^ Br (p,s)^l n (fx,f 2 ,f 3 ,xx ) x 2 ,x 3 ). (57) 

Here ^>i n is the light-cone wave function of the valence 
|3g) Fock component of the projectile proton, and \& / 
is the final state wave function of the recoil \3q) system 
after radiation of the QQ. 

The effective dipole cross section has the following 
structure, 



r? [rf a x (r 12 , p) + if ax (f 13 , p)] 



+ ^ fabc^ 1 T« [7f 2 (?X2 , P) + if <J 2 (f 13 , p)] 

+ quark permutations , 



(58) 



where the upper index g$ (i = 1,2,3) indicates the active 
quarks. 

The following new notation is introduced here, 

(Ti(r,p) = cr(f + ap) - a[f) , 
a 2 (f,p) = 2a[f—(l — a)p\ — a(f+ap)—a(f), 
rx2 = fx - f 2 , fx3 =fx-f 3 . (59) 

The differential cross section corresponding to the am- 
plitude Eq. (|57|) reads, 



daBripP ~* QQXp) 

dt' dzx dz 2 



J f 



One can employ completeness for the final state |3q) 
system,, 

^2 > r 2 , r 3 , xx , x 2 , x 3 )$j (r{ , f 2 , r 3 ' , x[ ,x' 2 ,x' s ) 
f 

3 

= \{5{n-f[)K^-^) ■ (61) 



Applying this condition to (f6"0")) and averaging over color 
indexes we get, 



d(JBr{pp — ► QQXp) 9 



dt' dzi (iz2 



256vr 



d 2 pd 2 s JJ^ d 2 ri 



2 i2 

x Tr Q ( s s P r(^ £ ri> ^2, F3) S^. f (p, s, fi, f 2 , r 3 )^ . 

(62) 

Here the trace is performed only over color indexes of the 
QQ pair, while the brackets indicate averaging over the 
colorless \3q) state in the proton. This trace results in, 



(p, s, fx, f 2 , f 3 ) S^ p r f (p, s, fi, r 2 , r 3 )) 



109 l I 

°i ( f i2 . p) + -^ a \{fvi,p) - ^^1(^12,^)^2(^12,/?) 

^ 0"! {fx2 , A) 0"! (fl3 . P) + g CT2 (fl2 , P)0- 2 [f 13 , p) 
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o - ! (^12, P)cr 2 {f l3 ,p) 



+ {fi2 t=± fi 3 } + quark permutations . (63) 
Further integrations in Eq. (|62p are straightforward. 



B. Leading twist production mechanism 

In this case the main contribution emerges from the 
diffractive interaction of a separate valence quark, al- 
though the interaction of spectators should be included 
as well. The full set of graphs having the leading twist 
behavior, 1/ttIq, in the cross section is shown in Fig. [7]). 
In all these graphs two gluons are attached to the target 




FIG. 7: Diffractive production graphs providing 1/ttiq con- 
tribution to the cross section. The projectile valence quarks 
should be permuted. 

and beam protons. The former is obvious, since we want 
to have a large rapidity gap and also that the recoil target 
proton remains intact. As for two gluons attached to the 
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beam, this condition is due to the leading twist behavior 
we would like to have. Indeed, in the set of graph where 
only one gluon is attached to the projectile, the second 
i-channel gluon is probing the small sizes r ~ 1/toq in 
the dissociation process G — * QQ. Thus one gains an 
additional suppression 1/toq. i.e. a higher twist behav- 
ior. Notice, however, that the pairs of gluons attached to 
the target and the beam are not symmetric, the former 
are in a colorless state, while in the latter all colors are 
summed up. 

The cross section corresponding to the graphs depicted 
in Fig. [7] has a form similar to (|62[) . 



— 3 

daprjpp -> QQXp) 9 f 2 2 tt 2 
dFd^ 2 = 256^ J d Pd S H d ndXl 

2 i2 

x ®in{ri,r 2 ,r 3 ,x l ,X2,x 3 ) $p r (p, s)\ 

x Tr Q (p p p r {p, s, ri, r 2 , r 3 ) Sp P r 1 (p, s, ri, r 2 , r 3 n . 

(64) 

In this case the effective dipole cross section has the form, 



+ 7* (rf ) [rf a 5 (f 12 , p, 5) + 7f a 5 (f 13 , p, 5) 

+ 7* (t^t?) [rf a 6 (f 12 , p, s) + rf ^(rW, 5) 
+ quark permutations . 



(65) 



where 



CT3 (p,s) 


= ^f(pv* 


(p - /3S) - <r[(l - a)p] 




+ cr[(l-a)p- 


f- (1 - /9)S] } 


0-4 (Pi S) 


= ^ftP + (l 






- o"[(l - ")P- 


-M} ; 




= a[ru - (1 - 


a)p] 




- o-[rij - (1 - 


a)p + l3s\ ; 


0"6 (rLj, p, S) 


= <r|f 14 - (1 - 


a)p - (1 - (3)s\ 




- a[r H -(l- 


")p1 , 


it 


= 2, 3 . 


(66) 



has the form, 

Tr Q ( E ?v (P) s , n,r 2 ,r 3 ) £f? r f (p, s , fi , f 2 , f 3 i 
19 2 2 

= y 0"3(Pi + 3 °"4(P) S) - - CT 3 (p, a) cr 4 (p, s) 
+ \ ai (P' [°" s ' P> + CT 5 (^13 , P, S)] 
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<T 4 (pi s) [cr 6 (fi 2 , p, S) + <T 6 (f 13 , p, s)] 



03 (p, s) [cr 5 (ri 2 , p, s) + 05(7=13, p, a)] 



+ gC r 3(p,s)[c r 6 (r ; *i2,p,s) + cr 6 (ri 3 ,p,s)] 

+ g [^5 (r i2 , p, ^ cr 6 (ri 2 , p, ^ + CT 5 (ri 3 ,p,s)cr 6 (fi 3 ,p,s)] 
5 

+ [^I(ri2,p,^ + ^(r 13 ,p,s) + (rl(ri2,p,s) + 
2 

+ °6 (^13 > P> S)] + g CT 5 ( 7 "l2 ; P, S) CT 5 (fl 3 , P, s) 

- ^cr 6 (fi 2 ,p, s)cr 6 (ri 3 ,p, s) - - [cr 5 (ri 2 ,p, s)<T 6 (ri 3 ,p, s) 
3 6 



(67) 



The dipole cross sections Ofc (k = 3,4,5,6) vanish lin- 
early in s as s — > 0. Since the mean value of s, con- 
trolled by the distribution function typ r (p, s), is small, 
(s 2 ) ~ 1/mg, one can make an expansion, 



er(r + si |2 ) — er(r) f» s*i,2 • V <j(r) 



(68) 



where si = 0s, s 2 = — (1 — /3)s, and r is one of the 
distances p, r^, etc. 



VI. SATURATION OF UNITARITY AND 
BREAKDOWN OF QCD FACTORIZATION 

Factorization assumes that the hard interaction of 
partons and subsequent hadronization proceed indepen- 
dently of the soft spectator partons in the beam and tar- 
get. This cannot be true for diffraction associated with a 
large rapidity gap (LRG). Indeed, the short range hard 
interaction of partons guarantees an overlap of the collid- 
ing hadrons, large impact parameters do not contribute. 
It is known from data that for such near central colli- 
sions unitarity is almost saturated [l6|, i.e. the chance 
for colliding hadrons to escape without soft inelastic in- 
teractions which terminate the LRG is very small. 

Here we rely on a simple eikonal model [6j, [28J, |35| ■ The 
absorptive corrections to the hard diffractive amplitude 
lead to a suppression factor, 



^(6)^^/(6) [l-ImfJQ>)] 



(69) 



The trace of the product of the effective cross sections 



where f%f(b) is the partial elastic amplitude. We as- 
sume a Gaussian shape for the elastic and diffraction 
amplitudes. After squaring the amplitude Eq. (|6l)|) and 
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integrating over impact parameter we arrive at the fol- 
lowing suppression factor for the diffractive cross section 
Eq. §1, 



K 



1 

1 



(s) 



7T B sd (s) + 2B p J(s) 



t' J ' J 



(4tt)2 B™(s) [B sd (s)+B™(s)] 



(70) 



Here the elastic slope depends on energy as B^f(s) = 
B ei+ 2a TP InO/so) with B% = 7.5 GeV -2 , s = lGeV 2 . 
The slope of single-diffractive hard cross section can be 
estimated as, B s d(s) « (r 2 /l )/3 + 2 ajp, where the proton 
mean charge radius squared is (r^ h ) — 0.8 fm 2 . 

A more accurate estimate needs a detailed information 
about the transverse structure of different Fock states. 
Such information is very much model dependent. Upon 
reaching the unitarity limit (Froissart bound) the fraction 
of diffractive events is expected to vanish as l/ln(s) [J]. 
How soon it may happen depends on specific models. For 
instance if gluons form dense spots inside hadrons, those 
spots can approach the unitarity limit much faster than 
the whole hadron-hadron scattering amplitude. Such 
black spots will suppress diffractive gluonic reactions 
(heavy flavors, triple-Pomeron term, etc.) much more 
than is suggested by Eq. (f?T)|) . Therefore, the predicted 
energy dependence of the survival probability Eq. ([70]) 
might be quite wrong and the diffractive cross section at 
the LHC energy may be overestimated. 

Gribov corrections [l2| were introduced into the sur- 
vival probability by means of a multi-channel treatment 
of the unitarity corrections in [3(| H3, [M, [39[ . The result- 
ing suppression factor is close to the eikonal one Eq. ([TO")) , 
differing by less than 10%. As an effective description, 
one can also incorporate the unitarity effects into the 
renormalized Pomeron flux [40] | . This description has the 
advantage of simplicity. 



VII. NUMERICAL RESULTS 

Now we are in a position to perform the integrations 
in Eqs. (|62| and (I64[) . Integrating the proton light-cone 
wave function squared, |\&i n (rj, %i)\ 2 , over all variables, 
except one of the x% , one gets the valence quark distribu- 
tion function, ^[2u v (x) + d v (x)]. Assuming a Gaussian 
dependence on r*j and factorized dependence on both 
and Xi , one can perform the full integration in (|62p , (|64|) 
and single out the contributions of the different mech- 
anisms. We rely on the phenomenological dipole cross 
section which has a saturated shape, 



a(r, x) = <T 



1 



(71) 



where the parameters have been fitted to HERA data 
for the proton structure function at small x in Ref. [4l| . 



First, we found that the bremsstrahlung mechanism, 
although leading twist, is very much suppressed. It con- 
tributes only a few percent of the production mechanism 
at the energy of RHIC, and order of magnitude less than 
that at the energy of LHC. Therefore, we can safely ne- 
glect the bremsstrahlung term. 

Second, we can disentangle between the contribution 
coming from diffractive excitation of individual valence 
quarks corresponding to the upper line of graphs in 
Fig. [7]) , and interference terms indicated by the bottom 
line of graphs in Fig. [71 The latter is controlled by the 
parameter Z = {r 2 h } p / Rq(x) where the numerator is the 
proton mean charge radius squared, (r^ h ) p ~ 0.8 fm 2 , 
and the denominator is defined in (|7ip . The interfer- 
ence terms vanish like 1/Z at large Z. This is because 
the diffraction amplitude is proportional to the difference 
between the cross sections of fluctuations of different size 
However, at r > Rq the cross section levels off and 
the interference diffractive amplitudes vanish. This gen- 
eral feature of diffraction is realized in Eq. (j68f . At high 
energies R (x) = 0.4 fm (i/xo) 0144 with x n = 3 x 10~ 4 
[4ll ] . Thus, for charm production i? ~ 0.35 fm (Z = 7) 
at the energy of RHIC, and R Q = 0.07 fm (Z = 160) at 
the energy of LHC. Numerical calculations using Eq. (j64|) 
confirm the smallness of interference terms, which pro- 
vide only about 1% of the cross section. 

Thus, with good accuracy we can neglect the inter- 
ference terms in the production mechanism. This step 
considerably simplifies the further calculations. Indeed, 
since the diffractive pp cross section comes out as a sum 
of diffractive excitations of the proton constituents, we 
can add sea quarks and gluons as well, i.e. make a re- 
placement, 



J2q{x) +q(x) + ^g(x) 
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(72) 



Here x = Attiq/xs. 



We remind the reader that diffractive excitation of a 
gluon should be calculated differently from that of a 
quark, as described in Sect. IIVI 

Next, we should specify the QCD couplings a s in (|45[) . 
(|5"Tj) . One of them, a s (l/s) corresponds to the hard scale 
of the reaction, s ~ 1/2toq. We use one loop approxima- 
tion with three, four and five flavors for the production 
of charm, beauty and top respectively. The coupling in 
(|4"S"]) and (f5"Tj) , a s (l /rho), should be taken on a soft scale, 
(1/p). As we have just explained, very large distances p 
are suppressed, since the saturated dipole cross section 
levels off and is independent of p. Therefore the typi- 
cal scale for this coupling is controlled by the saturation 
scale Q s = 2/Rq{x). Since it partially cover rather low 
values of Q s , the problem of infra-red behavior of a s may 
become an issue. We free ze t he coupling at the critical 
value gi| a s = 3^/4(1-^2/3) (see discussion in [HHl]. 

The results for the cross section of diffractive produc- 
tion of charm, beauty and top, p + p — > QQX + p, are 
plotted as function of energy in Fig. [5J To be compared 
with available data (see next section) charm diffractive 
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FIG. 8: Cross section of diffractive production of heavy flavors 
as function of energy. The experimental points are the results 
of the E690 and CDF [3 experiments. 



cross section is integrated over xf > 0.85, and beauty 
over xf > 0.9 (same for top). All the cross sections 
steadily rise with energy. The cross sections of charm 
and beauty production differ by about an order of mag- 
nitude what confirms the expected leading twist behavior 
l/m 2 Q . 

We also calculated the x\ distribution of a diffractively 
produced charm quark by integrating over all other vari- 
ables. X\ = /pp is the ratio of plus components of the 
produced c-quark and the incoming proton. The results 
are shown in Fig. [5] at RHIC and LHC energies. 

Notice that to be compared with data (unavailable so 
far) for production of charmed mesons, this result has 
to be corrected for the fragmentation c — ► D which is 
poorly known. The resulting behavior at X\ — > 1 should 
obey the end-point behavior dictated by Regge. There- 
fore we expect it to be less steep than what is plotted in 
Fig. [5J One may wonder: a convolution with the frag- 
mentation function c — > D may only result in a steeper 
fall off at a; i — > 1, how can it become less steep? The 
answer is: the convolution procedure is incorrect, QCD 
factorization badly fails at x\ — > 1. The usual fragmen- 
tation function measured, say, in e + e~ annihilation, cor- 
responds to a fast c-quark producing a jet and picking up 
a slow light quark from vacuum to form a D-meson. In 
hadronic collisions at large x\ hadronization occurs dif- 
ferently: a fast projectile light quark picks up a slow o 
quarks produced perturbatively. Correspondingly, in the 
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FIG. 9: The cross section of diffractive excitation of a proton 
with charm production as function of fraction x\ of the initial 
proton momentum carried by the charm quark. The cross 
section is shown at the energies of RHIC and LHC. 



case of diffractive production of a heavy flavored baryon 
a leading projectile diquark can pick up the heavy quark. 

Notice also that X\ has a bottom bound imposed by the 
kinematics of diffraction, x\ > 4mg/(l — xp)s, where Xf 
is the Feynman variable of the recoil proton in pp — > Xp. 
In order to comply with available data (see next section) 
we integrate over xf > 0.85 for charm (also top), and 
xf > 0.9 for beauty. 

Our results for transverse momentum distribution of 
diffractively produced quarks are presented in Figs. (fT0|) - 
(fP2|) for different heavy flavors and energies. 

There px distributions hardly correlate with x\ of the 
heavy quark, what is quite different from the usual sea- 
gull effect. We remind, however, that this is not the 
usual factorization based hadronization. In this case a 
fast projectile quark-spectator picks up a slow heavy fla- 
vor. Therefore, the transverse momentum of the pro- 
duced heavy flavored meson is mainly controlled by the 
transverse momentum of the light spectator. 

To conclude this section, we should comment on the ac- 
curacy of performed calculations. The main uncertainty 
seems to be related to the absorptive (unitarity) correc- 
tions. Comparing different models, the difference is not 
dramatic, of the order of 10%, with a probability fac- 
tor K = 0.14 at the Tevatron energy. However, all those 
models may miss the specific dynamics of interaction dis- 
cussed in Sect. IVII and overestimate diffraction at the 
LHC energy by much more than 10%. The next theoret- 
ical uncertainty is related to the choice of heavy quark 
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FIG. 12: The pr dependence of the cross section of diffractive 
excitation of top. 
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FIG. 11: The pr dependence of the cross section of diffractive 
excitation of beauty. 



masses. We used m c — 1.5 GeV, m& = 4.7 GeV and 
m t — 175 GeV. Our diffractive cross sections approxi- 
mately scale as 1/toq. The cross section also depends on 
the behavior of the QCD coupling in the infrared limit. 
We freeze the coupling at the critical value a s < 0.4. A 
larger value would lead to a corresponding increase of the 
cross section. 



VIII. DATA 
A. Diffractive production of charm 

Data for diffractive production of heavy flavors are 
scarce. Three experiments searched for open charm pro- 
duction in diffractive events. Diffractive production of 
A c was first detected in experiment [45[ at y/s — 63 GeV 
at ISR. The reported cross section u{pp — > A c Xp) = 
10 — 40/ib is quite above our calculations. However, 
one should take this experimental result with precau- 
tion, since it also considerably exceeds the results of a 
later measurement (see below). The same experiment 
overestimated the total cross section of inclusive charm 
production by more than an order of magnitude (see dis- 
cussion in [44j|). 

A stringent upper limit for diffractive charm produc- 
tion at y/s = 40 GeV was imposed by the E653 experi- 
ment at Fermilab [46| . which found the cross section to 
be smaller than 26 /ib for proton-silicon collisions. These 
measurements included events with the nucleus remain- 
ing intact, as well as fragmenting. The nuclear effects are 
controlled mostly by the survival probability of the large 
rapidity gap during propagation of projectile soft spec- 
tator partons through the nucleus. Therefore, only nu- 
clear periphery contributes and nuclear effects are rather 
small. According to calculations in [H, |47j, the cross 
sections on silicon and free proton are similar, so one can 
apply the measured upper bound for pp collisions as well. 
Our calculations are well below this bound. 

An accurate measurement of the cross section of 
diffractive production of D*-mesons was performed in the 
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E690 experiment at Fermilab at s/s = 40 GeV and xp > 
0.85 [44|. The result for the cross section of cc production 
based on the most accurate data for production of D*~ is 
<j{p + p^ ccX+p) = [0.61 ±0.12{stat) ±0.11{syst)} fib. 
This cross sections is integrated over xf > 0.85. Plotted 
in Fig. [8] it agrees well with our calculations. 



B. Diffractive production of beauty 

An upper limit for the cross section of diffractive 
beauty production at \/s — 630 GeV was established by 
the UAl experiment, with a large theoretical uncertainty: 
<j(pp — > bX +p/p) < [0.6 — 1.2] fib. The corresponding 
upper limit for the fraction of diffractive production rel- 
ative to inclusive production of beauty was found to be 
R bb < [3.1 - 6.2]%. 

The first observation of diffractive beauty produc- 
tion was done by the CDF collaboration [48[ at y/s — 
1800 GeV and xf > 0.9. Unfortunately the cross section 
suffers large theoretical uncertainties. The Monte Carlo 
codes used for the evaluation of the gap acceptance were 
based on factorization, which is broken for diffraction. 
Depending on the assumed distributions of quarks and 
gluons in the Pomeron, the fraction of diffractive events 
with x f > 0.9 ranges from R bb = [0.62 ± 0.19(stat]± 
0.16(syst)}% to [1.18 ± 0.36(stat) ± 0. 27 (syst)]% gg[. 
From these results one can probably conclude only that 
the ratio is of the order of 1% with the error of the or- 
der of 0.5%. We use this estimate for further comparison 
with our calculations. 

It also worth mentioning that the electrons from b- 
decay are detected with large transverse momenta pt > 
9.5 GeV |48| . As far as diffraction this is leading twist, 
and it should have a p-r-dependence similar to that in in- 
clusive production. Then the experimental p^-cut should 
not affect much the diffraction-to-inclusive ratio. 

To get the diffractive cross section one needs to know 
the inclusive cross section of beauty production which 
unfortunately has not been measured at this energy. 
One can rely on the result of the UAl experiment at 
■/s = 630 GeV, which is a(pp -> b + X) = [19.3 ± 
7(exp) ± 9(theor)] fib, and also on data at lower ener- 
gies [50L IBll HH and a theoretical extrapolation. Two 
theoretical approaches explain reasonably well data for 
inclusive production of beauty, the NLO parton model 
[53[ ] and the light-cone dipole description (3l|. Although 
the absolute predictions differ by about a factor of 2, 
they predict the same energy dependence. With this en- 
ergy dependence we extrapolated the result of the UAl at 
xfs = 630 GeV to the Tevatron energy y/s = 1800 GeV, 
and found <r(pp — > bbX) = [103.5 ± 86] fib (we added 
linearly the experimental and theoretical uncertainties). 

This result is quite uncertain because of the large the- 
oretical error and of the long energy interval for ex- 
trapolation. Another possibility to make an estimate is 
to use CDF data [54[ for inclusive 6-quark production 
at \/s = 1960 GeV in the rapidity interval \y\ < 0.6, 



a = [17.6 ± 0A(stat)~^2'^(syst)] fib, and to extrapolate it 
to other rapidities. A rough estimate would be to assume 
the same production rate for the whole rapidity inter- 
val Ay = 2 \a.(^fs / rubr) ■ We take the mean transverse 
mass m1 T — m b + (p^) ~ 2m\. The l as t small correc- 
tion is to scale these numbers according to the theoreti- 
cal energy dependence [31] . <r(\/s = 1960GeU)/cr(v / s = 
1800GeU) = 0.95. Eventually, we arrive at the estimate, 
a(pp ^b + X) = [139.7 ± 3(stat)±H(syst)\ fib. 

Alternatively, one can rely on theoretical predictions 
for the inclusive cross section done within the dipole ap- 
proach. Usually such calculations are pretty accurate, 
since are based on the phenomenological dipole cross sec- 
tion fitted to DIS data, and include all higher order cor- 
rections and higher twists. In particular, it describes 
quite well the available data for inclusive production of 
charm. The predicted inclusive cross section of inclu- 
sive beauty production at \/s = 1800 GeV is 140 fib, 
which agrees well with the above extrapolated experi- 
mental cross sections. 

Relying on this theoretically predicted inclusive cross 
section and experimentally measured [48| fraction R bb 
of diffractively produced beauty, we finally arrive at the 
diffractive cross section, <J s d{pp ~ ► bXp) — [1.4 ± 0.7] fib. 
This value is plotted in Fig. [8] in comparison with our 
calculations. Still, one should remember that both the 
value and error are subject to considerable uncertainties, 
in particular theoretical ones. 

IX. SUMMARY 

The main results of this paper can be summarized as 
follows. 

• Novel leading twist mechanisms of diffractive ex- 
citation of heavy flavors in hadronic collision are 
proposed and calculated. Factorization leading to 
higher twist diffraction is badly broken. 

• Two mechanisms of heavy flavor production are 
identified (see Fig.[T]). One, called bremsstrahlung, 
is similar to the Drell-Yan mechanism of radiation 
of a heavy dilepton, but includes also interaction of 
the radiated virtual gluon with the target. Another 
mechanism, called production, involves also inter- 
action of the heavy quarks with the target. Diffrac- 
tion excitation of QQ in a separate parton by the 
bremsstrahlung mechanism is higher twist and can 
be neglected. The leading twist excitation of pro- 
jectile quarks or gluons is possible only via the pro- 
ductive mechanism. 

• The presence of spectator partons in the projec- 
tile hadron opens new possibilities of interactions, 
and the bremsstrahlung mechanism becomes lead- 
ing twist as well. Quantitatively, however, it is still 
a small part of the cross section. The dominant 
contribution comes from diffractive excitation of a 
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separate projectile quark or gluon via the produc- 
tion mechanism. 



Available data for diffractive production of charm 
and beauty agree with our calculations. The lead- 
ing twist dependence on the quark mass is also con- 
firmed. 



• Our results for heavy flavors allow straightforward 
application to high-p^ jets. The same leading 
twist production mechanism explains the observed 
independence of hard scale of the diffraction-to- 
inclusive ratio for di-jet production [55[ . Numerical 
calculations and comparison with data will be done 
elsewhere. 
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Appendix A. CLASSIFICATION OF THE 
AMPLITUDES 



Ml 



Mi 



Mi = 



M 



Mi 



E^(-k)^)r^r^r^ f a (q) 

A . 

[M 2 (l - a) + a 2 m 2 q + k 2 ][m 2 Q + k 2 } ' 

Er^r^K 1 f a (q) 

A . 

[M 2 (l - a) + a 2 m 2 q + n 2 ][m 2 Q + k 2 } ' 

E r£ (p 2 ) r 2 A (k) if abc 7* n j a (q) 

_a . 

\p\ + a 2 m 2 ] [niQ + n 2 ] 

E rK&) r*(* - m rf r b Q T Q fM 



(A.l) 
(A.2) 
(A.3) 



[p 2 + a 2 m 2 q ] [m 2 Q + (k - (3q) 2 + (3(1 - (3)Q 2 ] ' 

(A.4) 

1 

p 2 + a 2 m 2 

ErH-P2)T^K+(l-mr^rQr^ f a (q) 



m 2 Q + [k + (1 - + 13(1 - P)Q 2 



(A.5) 



Here, f a (q) is the amplitude for emission of a gluon with 
color index a and transverse momentum q by the target 
proton; 

TiW = X/ {(2-a)(e A ■ tt) + ia[a x e A ] ■ n 

+ ia 2 m q [a x e A ] ■ n| Xin, 
r*(FZ) = (j> ] [mi{? ■ e x ) + (1 - 2/3) • n)(e A • Ft) 



+ i\e\ x n 



(A.6) 



Notice that within each of the two groups of ampli- 
tudes, Mj , Mj, and Mj, Mj, the denominators have 
similar structure and these amplitudes may be combined 
producing light-cone distribution amplitudes in impact 
parameter representation. Only the amplitude M^ does 
not belong to any of these groups. Nevertheless, it can 
be splited into two parts. The first one comes out after 
multiplying Mj by factor, 



p\ + a 2 m\ 



M 2 + Q 2 M 2 (1 - a) + a 2 m 2 + p 2 



(A.7) 



Then the amplitude acquires the denominator of the first 
group, which we call bremsstrahlung. The rest of M-f , 
which gets the factor, 



A.l. Transverse polarization 



The five amplitudes corresponding to the Feynman 
graphs in Fig. [21 with transversely polarized gluons 
radiated by the projectile quark have the following form, 



M 1 



M 2 + Q 2 



i 2 Q + k 2 + (3(1 - (3)Q 2 ' 



(A. 



has the structure corresponding to the second group, 
which we call production mechanism. This explains the 
way in which we classify the amplitudes in Eqs. I©-©. 
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Using the relation, 

(A.9) 

we obtain for the bremsstrahlung and production ampli- 
tudes Eqs. ([8])- ([9]), respectively, 

M Br = {[®l r (T? + (I - a)q, k) - <t>l r (n - aq, k)] 



a 'b 'b 



- $ 



Br\ 



(A.10) 



where 



$l r (7T, K) 



£1^)1^) 

A 



[M 2 (l - a) + a 2 m 2 q + t? 2 ] (m 2 Q + k 2 ) 



(A.ll) 



Correspondingly, for the production mechanism, 
Ml r = { [*pt (P2 ,/e-/3g) - $L (ft , «)] Tf r^rf 
+ $J r (p2 ! K)-$J r (p 2! /t+(l-/3)g) 



x 7? r: 



(A.12) 



ErK^)r 2 A («) 

A 



where 

$£ r (p 2 ,/c) = 

(pi + a 2 m 2 ) m| + k 2 + 0(1 - 0)Q 2 

(A.13) 

Notice that the amplitudes Eqs. (|A.10p , (|A.12|) vanish 
in the forward direction, at q — ► 

Now we can convert the distribution amplitude 
Eq. (|A. 131 to impact parameter representation, 



Br(Pr 



(A.14) 



7a(&) = ^ / d^e^ f a (q) • 



(A.15) 



Then the amplitudes Eqs. (|A.10|) . (|A. 12|) take the form, 

x { [la (& - (1 - - 7a(& + a/5)] r«7f 7? 
+ [7a(6)-7-(& + «^)] t?t?t?}; (A.16) 



A. 2. Longitudinal polarization 

The five amplitudes corresponding to the graphs in 
Fig. [2] 



4r> 6 % Q /„(£)(x|xin)(0 t <?-n 

1 



'6 '6 VA/ , 

1 -a 



_M 2 (l-a) + a 2 m 2 + K 2 M 2 
1 - a 1 



Mi- = 



Ml 



M 5 L 



M 2 (l-a) + a 2 m 2 + n 2 M 2 
Aifabc Tf rf / a (<f) ( X J Xvn) (0 f a • n , 
1 1 
^ ~ M 2 _ 

4r t 'V^ /a® (xj Xin) (0 f <? • n0) 

(3(1-0) l_ 

m 2 Q + (n- 0q) 2 + Q 2 0(1 - 0) Q 2 

ir^r^f a (q)( X \xin)(<P^-n^ 

0(1-0) 



1 



l 2 Q + (R-(l-0)q)+Q 2 0(1-0) Q 2 

(A.18) 



Summing up these amplitudes, the terms proportional 
to 1/Q 2 and 1/M 2 cancel. The rest, the amplitudes 
(i = 1, 2, 4, 5) without these terms, which we denote M^, 
can be grouped in a way to create light-cone distribution 
amplitudes. To reach this goal we introduce an additional 
amplitude which is identical to zero, 



where 



Mi = M 3 L 1 - M 3 L 2 = 



= 4i/ o6c 7f f a (q){x)^n)(<P" ff- 

1 - a 



(A.19) 



M 2 (l-a) + m 2 +p 2 ' 



M£ >2 = Aif abc T^T b fa(q)(x)Xr n )(4>^-n4>) 
0(1 - 0) 



m 2 + K 2 + Q 2 0(1 - 0) 



(A.20) 



d 2 bd 2 P d 2 s e m+i *P +iH - g $l r (p, s) 
{ [la (b- (1 - a)p + 0s) - Ja (b - (1 - a)p) 

la (b-(l- a)p*j 
7a (b-(l-a)p-(l-0)s)] t^t^t?}. (A.17) 



x T b T b T a 



Now we are in the position to group the longitudinal 
amplitudes in a similar form to the transverse ones, get- 
ting the light-cone distribution amplitudes corresponding 
to the bremsstrahlung and production mechanisms, 



M| r = Mi 



M 2 +M- 



3,1 ) 



M^ r = Mi + Mr, — M 3 L 2 . 



(A.21) 
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These amplitudes have the form, More relations for r-matrices, 



M| r = 



{ [* L Br (7? + (1 - a)<f, S) - $ L Br (vr - ag, «)] 7* r b 91 t£ 
{ [$| r (vf, k) - $^(7? + (1 - a)£ k)] T«W«W?}f a (q); 



Trr« = 0; Tr [r a r b ] = £<U; 
Tr [r Q T b r c ] = - /i abc , ft, abc = d abc + if a bc', 



(A.22) 11 L a c J - 4 



Tr [T a T b T c T d ] = — 6 ab S cd + - h a be h c de, (B.2) 



[$| r (p 2) *) - $L(P2,K + (1 - /?)$)] « r b Q }/ a (q), 

(A.23) 



where the longitudinal distribution amplitudes <& read, 

r 4(l-a)(xtxm)(0 t CT-n0) r Q r b r a = ~n; 

BrK ' M 2 (l - a) + a 2 m 2 q + ir 2 1 1 



t q Ti, r c r a = - S bc ~ -T b T c ; (B.3) 
4/3(1 -/?)(x^Xin)(^t?-^) 4 6 

+ k 2 + Q 2 /?(l - /3) ' 

(A.24) 



faklfbkl — 3<5 ab ; 

Appendix B. USEFUL ALGEBRAS 5 

dakldbkl — g <5ab; 



Average of a product of arbitrary functions , , , 3 

A(r), B(t), C(t) of Gell-Mann matrices over the Jiakhbltld - 



proton wave function has the general form, 3 

Uiakjkbljlci ~ 



2 ^aoc j 



(A(r*)B(r«)C(r»)) N 5 

1 ^ ^iak^kbl J lei — ^ Jabc-) 

-|TrA(r)Tri?(T)TrC(r) J 

Tr [A(r) B( r ) C(t)] + Tr [A(r) C(r) B(r)] d iak d kbl d lci = - d abc . (B.4) 
Tr A(r) Tr [B(t) C(t)] - TrB(r) Tr [A(r) C(r)] 
TrC(r)Tr[A(r)i?(r)]} (B.l) 
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